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H Abstract 

We investigate the inclusive production of the Higgs bosons of the minimal supersym- 
metric model with a heavy quark pair, (tt and bb), in e + e _ collisions at high energies, 
y / s = 500 GeV and y / s = l TeV. The 0{a s ) QCD radiative corrections are included and 
the dependence of the production rates on the parameters of the MSSM is explored. The 
associated production of a Higgs boson with a bb pair receives large resonant contributions 
and can have a significant rate. 
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Figure 1: Mass of the lightest neutral Higgs boson as a function of the pseudoscalar mass 
Ma- The squarks are assumed to have a common mass, Ms = 500 GeV and the scalar mixing 
parameters are set to zero. 

1 Introduction 

The search for the Higgs boson is an important objective of present and future colliders. In the 
simplest version of the Standard Model, the Higgs boson is necessary in order to understand 
fermion and gauge boson masses. Once the Higgs boson has been discovered, it is crucial 
to measure its couplings to fermions and gauge bosons. The couplings to the gauge bosons 
can be measured through the associated production processes, e + e~ — > Zip, qq' —> W ± (p, and 
qq — ► Z(j), and through vector boson fusion, W + W~ —>■ <j) and ZZ — > </>, (where is a generic 
Higgs boson). The couplings of the Higgs boson to fermions are more difficult to measure, 
however. 

In the Standard Model, the fermion-Higgs couplings are completely determined in terms of 
the fermion masses, 

M Q m 
9QQh.su — > \ l ) 

where v = (y^Gp)^ 1 ^ 2 and so the top quark-Higgs boson Yukawa coupling is large while the 
bottom quark-Higgs boson Yukawa coupling is much smaller. In extensions of the Standard 
Model, however, the Yukawa couplings can be significantly different. They are no longer strictly 
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proportional to the fermion masses, but depend on the parameters of the model. The mini- 
mal supersymmetric model (MSSM) provides a useful benchmark for comparison of the QQ<p 
couplings with those of the Standard Model. 

In fact, the process e + e~ — > ti<j) provides a direct mechanism for measuring the top quark- 
Higgs Yukawa coupling, and for discriminating the Standar Model type Yukawa coupling from 
others. |],| This process proceeds via 7 and Z exchange, with the Higgs boson bremsstrahlung 
from the heavy quark. There is also a contribution from the Higgs boson coupling to the 
exchanged Z, but this is always less than a few percent and so does not upset the interpretation 
of this process as a measurement of the top quark-Higgs Yukawa couplings. 

The dominant production mechanism for a Higgs boson at a hadron machine is gluon fusion, 
gg — > 0, which proceeds by a loop diagram containing heavy fermions and so is also sensitive to 
the fermion-Higgs Yukawa couplings. For heavy fermions, however, the gluon fusion process is 
independent of the fermion mass and so this mechanism counts the number of heavy fermions 
which couple to the Higgs boson. The process gg — > <p thus does not provide an unambiguous 
determination of the coupling of the Higgs boson to heavy fermions since it is sensitive to the 
unknown number of heavy fermions. In addition there are uncertainties due to the fact that 
QCD corrections to the gg — > h process are large. 

In a previous paper ||, we computed the 0(a s ) corrections to the process e + e~ — > tth SM . 
At y/s = 500 GeV, they are large and positive, while at y/s = 1 TeV, the QCD corrections 
are small. In this paper, we apply the 0(a 3 ) corrections |3],(4] to the production of the Higgs 
bosons of the MSSM and explore the dependence of the production rates on the parameters of 
the MSSM. 

It is more difficult to extract the bottom quark-Higgs Yukawa coupling from a measurement 
of e + e~ — > bbhsMi since the coupling itself is tiny and the Z contribution is relevant. However in 
the MSSM, for certain values of tan /3, the ratio of the neutral Higgs vacuum expectation values, 
the bb<f) Yukawa coupling receives significant enhancements and so the process e + e~ — > bbcf) may 
be larger than in the Standard Model. Moreover, there are 5 Higgs bosons in the MSSM, so 
that additional processes not present in the Standard Model may be useful to pin down the 
fermion-Higgs boson Yukawa couplings. In particular, there are new resonance contributions, 
for example e + e _ — ► A°hi,A° — > bb, which are not present in the Standard Model. The effects 
of these new contributions can be significant and are the focus of this paper. 

The process e + e~ — > bb<j) is not only sensitive to the b-quark Higgs Yukawa coupling but also 
has an interesting dependence on the other parameters of the supersymmetric model. Sections 
5 and 6 explore this dependence in detail. 

In Section ||], we provide an overview of the elements of the MSSM relevant for our study. 
The following sections describe our results for the associated production of the scalars and 
pseudoscalar of the MSSM with it and bb pairs. Section [7| contains some conclusions. 
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2 SUSY Background 



The minimal supersymmetric version of the standard model (MSSM) is a theoretically appealing 
extension of the Standard Model and provides a useful framework for comparison with the 
Standard Model. In the MSSM, the Higgs sector contains two SU{2)l doublets and so there 
are five physical Higgs bosons; two neutral scalars, h° and H°, (which we will often denote by 
h®), a pseudo-scalar A , and a pair of charged scalars, H ± . || We are therefore led to consider 
the processes, 



e + e~ -> QQh° , 
e + e~ -> QQH° , 
e + e~ -> QQA° , 
e + e" -> M + , , 

(2) 

where Q is either t or 6. These processes have been previously considered in Ref. [0 and we 
extend their analysis. The rates are small, but the signatures are distinctive.^] 

At lowest order, all masses and couplings in the Higgs sector can be computed in terms of 
two parameters, which are typically taken to be Ma, the mass of the pseudoscalar Higgs boson, 
and tan (3, the ratio of vacuum expectation values. At higher order, there is a dependence on 
the scalar masses and mixing parameters. Having picked Ma, tan/?, and the scalar masses and 
mixing parameters, the neutral Higgs boson masses are predicted [H-f|. 

The MSSM has the feature that there is an upper limit on the mass of the lightest neutral 
Higgs boson. Including the radiative corrections to leading logarithm, we have, 

Ml < M 2 Z cos 2 2(3 + 5 , (3) 

with 



3G F M* ( M 



V^TT 2 b V M, 



5="-^- log l + £f , (4) 



and Ms is a common squark mass. In our numerical studies the renormalization group improved 
values for the Higgs boson masses, along with the NLO corrections are included, using the 
program HDECAY. [[| Fig. [I] shows the mass of the lightest neutral Higgs boson as a function 
of Ma- There is clearly a maximum value around Mh = 110 GeV. Including mixing effects in 
the squark and Higgs sectors raises this limit to around 130 GeV. Such a light Higgs boson 

1 The processes e + e~ — > tbH + , tt>H~ are interesting in the small mass region where i? ± cannot decay to t 
nor t decay to i/ ± [0]. Outside this mass region, other processes are more useful for determining the Yukawa 
couplings. (|] We will not consider these processes further. 
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Table 1: Higgs boson couplings to fermions in units of gQQh 8M = ~Mq/v 
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may be observable in association with a heavy quark pair, providing additional impetus to our 
study. 

The LEP2 experiments have searched for the processes e + e~ — > and e + e~ — > h°A° and 



find the 95% confidence level limits, (10 



M A 
tan (3 



> 78 GeV 

> 2.1 , 



(5) 



for the case with no mixing in the squark sector .0 We will therefore consider tan (3 = 2.5, 5 and 
40 and Ma > 80 GeV0 as representative values in our study.Q The corresponding values for the 
light Higgs mass, Mh, can be found from Fig. [I]. 

An important feature of SUSY models is that the fermion-Higgs couplings are no longer 
strictly proportional to mass as they are in the Standard Model. It is convenient to write the 
couplings for the neutral Higgs bosons to the fermions in terms of the Standard Model Higgs 
boson couplings. In particular, we write the couplings of the Higgs bosons to the quarks as 



CQQhQiQih + CQQ H QiQiH° + CQQAQii^QiA^ 



(6) 



where Cqq^ is 1 for a Standard Model Higgs boson, (Cqqh = Cqqa = 0), while the corresponding 
MSSM couplings are given in Table 1. The fermion-Higgs boson Yukawa couplings in the MSSM 
are then, 



QQQhi = -(V2G F ) 1/2 M Q C QQhi , 

9qqa = -(V2G F ) 1/2 M Q C QQA . (7) 

For small Ma, the couplings of the neutral Higgs boson to fermions can be significantly different 
from the Standard Model couplings; the 6-quark coupling to the lighter of the neutral Higgs 

2 There is also a small allowed region with tan (3< 0.8. Jl0| 

3 Results from precision electroweak measurements suggest that the true bound on Ma may be slighter than 
this. fuj. 

4 Of course, for a given value of tan/3, the lower bound on Ma may be greater than that of Eq. ||. 



4 




h°,H°A° 



+ 




Figure 2: Feynman diagrams contributing to the lowest order processes, e + e 
the MSSM. 
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Figure 3: Feynman diagrams contributing to e + e 
Standard Model. 



QQh®,A° which are not present in the 



bosons, h°, becomes enhanced, while the t-quark coupling to h° is suppressed for large tan/3. 
On the other hand, the coupling of the heavier neutral Higgs boson, H°, to the b is enhanced 
at small tan/3 and the coupling to t enhanced at large tan/3 for light Ma- 

The angle a is related to the mixing in the neutral Higgs boson mass matrix. In the leading 
logarithmic approximation, 

_ (M% + M|) sin 2/3 
tan ~ (Mi-Mf) COS 2/? + 5/sin 2 /? ' [ } 

where — 7r/2<a<0. In our numerical studies we include the renormalization group improved 
values for a using the program HDECAY. || 
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The rate for e + e 



ah®, h® = h°, H°, receives the contributions from 7 and Z exchange shown 
in Fig. ^| as in the Standard Model. In the MSSM, there is an additional contribution not 



present in the Standard Model from the sub-process e + e 



A hi, A 



tt, shown in Fig. |[ 
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The rate for the associated production of a neutral SUSY Higgs boson with a heavy quark 
can be written as J2],[| 



da{e+e- - QQh ,) N a 



dxh (47r) 



v\ + a 2 



2 2 2Q e Q Q v e v Q (^ 2 + Qe)K + ag) 
V e Vg + 1-Mf/s + (1-Mf/s) 2 

G 7 



QeQQVeVQ \ . . 

+ 1-M1/, 6 / ' (9) 

where cr = 47ra 2 /3s, a is the QED fine structure constant, iV c = 3 is the number of colors, 
Xh = 2E h / \/s with the energy of the Higgs boson and Q;, and a, (i = e, Q) denote the 
electromagnetic and weak couplings of the electron and of the heavy quark respectively, 

Vi = 2/ 3L - jQiSw ^ a , = 2/ 3L ^ QqX 

with I\ L = ±1/2 being the weak isospin of the left-handed fermions and s 2 ^ = 1 — c 2 ^ = 0.23. 
The contribution from the ZA ^ coupling shown in Fig. |^ and its interference with the other 
contributions is given by the factor G 7 . For completeness we list the coefficients Gi(xh) in 
Appendix 

In the case of tth° and ttH° production, the dominant contribution at = 500 GeV is 
photon exchange, as is the case also in the Standard Model. In Fig. f|, we show the ratio of the 
cross section computed including only the photon exchange to the total cross section for the 
case of e + e~ —>■ tih°. Very similar results hold for the e + e~ — ► tiH° process. Since the cross 
section is dominated by photon exchange, the contributions from e + e~ — > A°h® — > tth® and 
e + e~ — > Z — > tthl can safely be neglected at yfs = 500 GeV. 

The 0(a s ) radiative corrections can be summarized by a K factor multiplying the lowest 
order rate. The K factor can be defined as 

. (11) 

From Fig. [|, it is clear that the K factor computed in Ref. ||, where only the 7 exchange 
was considered, can be used at ^/s = 500 GeV. The K factor depends only on the mass of the 
produced Higgs particle, and not on the SUSY parameters in the limit in which only the photon 
exchange contribution is important and so the results for the Standard Model Higgs boson can 
be used to obtain the 0{pt s ) corrections. The only dependence on the renormalization scale \l 
is in the strong coupling constant, a, s [p). We will take fi = M t = 175 GeV everywhere. A useful 
numerical fit to the K factor of Ref. H], valid for -^5 = 500 GeV, is, 
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Figure 4: Dominance of the photon contribution to the process e + e — > tth at ^ = 500 GeV. 



The squarks are assumed to have a common mass, 
parameters are set to zero. 



Ms = 500 GeV and the scalar mixing 



K U=50oGeV~ l + « s (/i 2 ){ 



64 



where 



and 



{yf~s~M h 



100 GeV 



4M 4 2 



+ a + air + a 2 r' 



(12) 



(13) 



«o 
a i 
a 2 



-1.760 
-.004 , 
-.133 . 



(14) 



This fit (valid only at yfs = 500 GeV), is an excellent approximation to the exact K factor 
of Refs. UK]. Fig. ^| shows that the QCD corrections increase the rate significantly at yfs = 
500 GeV. At yfs ~ 500 GeV, much (but not all) of the K factor can be accounted for by the 
threshold contribution to tt production. The first term in Eq. ([12]) is the threshold contribution 
to the K factor [ffl, which is plotted as the curve -ft'thres/io/d i n Fig. ^, and the rising of the K 
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Figure 5: K factor for the process e + e — > tthgM at ^=500 GeV. This figure uses a s (M t ) = 
0.11164. K exact is the result of Ref. |§ and K t hreshoid includes only the it threshold contribution. 



factor signaling the Coulomb singularity as the threshold is approached, \fs ~ 2M t + M hsM , is 
clearly seen. 

We can now compute the rates for e + e~ — > till® production including the 0(a s ) corrections. 
In Figs. |] and [7|, we show the lowest order and the radiatively corrected rates. All contributions 
to the lowest order rate (Z exchange, A°h® production, etc) are included and the next-to-leading 
order result is computed by multiplying the lowest order rate by the K factor of Fig. || Since 
the dominant contribution is photon exchange, this should be an excellent approximation. The 
rate for ith° is enhanced at small tan f3, while the rate for itH° production is enhanced at large 
tan/3. The rate is never larger than a few fb, making this process extremely challenging to 
observe. 

Although the processes e + e~ — > ith° and e + e~ — > itH° have small rates, the signatures are 
distinctive since the final state will be predominantly W + W~bbbb and so may be observable 
with a small number of events. In Fig. we show the region of the Ma — tan (3 plane where 
the rates for either process are larger than 0.75 fb and we see that this includes much of the 
parameter space. We have everywhere set the mixing in the squark and Higgs sectors to zero. 
Our results are relatively insensitive to this assumption. 

At y/s = 1 TeV, the cross section is no longer dominated by photon exchange and the Z 
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Figure 6: Cross section for e + e~ — > production at yfs = 500 GeV at lowest order and at 
NLO for different values of tan/3. We take a s (M t ) =0.11164 and use the K factor of Ref. [§]. 
The squarks are taken to have a common mass, Ms = 500 GeV, and the mixing parameters are 
set to zero. 



exchange becomes important. The ratio of the cross section computed using only the photon 
exchange to the total cross section varies between 70 and 80 %. However, the contribution 
from the A°h® production of Fig. |3| is always less than 1 %. In Ref. [|J, the K factor for 
Standard Model Higgs production was computed including the effects of Z exchange. Using 
a s (M t ) = 0.11164, the result of Ref. || is that the K factor at yfs = 1 TeV is roughly a 
constant, K(lTeV) ~ 0.94. We will use this result in our plots. Since the radiation of a 
Higgs boson from the Z boson (the second diagram in Fig. is small, the K factor is again 
approximately independent of the SUSY parameters and depends only on the mass of the 
produced Higgs boson. At high energy where yfs is not close to the kinematic threshold, the 
dominant contribution to the QCD corrections is from the Higgs vertex correction computed 



in Ref. [12]. In this limit the K factor is roughly, 



^1_^£M^0.9, (15) 

7T 

in agreement with the result found from the complete calculation for a/s = 1 TeV. HU Figs. ||] 
and |T0| show the rates for e + e _ — ► tth® at y/s = 1 TeV. The cross section for either e + e~ — > tth° or 
e + e~ — > tiH° is greater than 2 fb throughout most of the Ma— tan (3 plane. For Ma > 120 GeV, 
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Figure 7: Cross section for e + e~ — > ttH° production at y/s = 500 GeV at lowest order and at 
NLO for different values of tan/3. We take a s (M t ) =0.11164 and use the K factor of Ref. 
The squarks are taken to have a common mass, Ms = 500 GeV and the mixing parameters are 
set to zero. 

the rate for e + e~ — > iih° is relatively insensitive to tan (3 and is always between 2.4 and 2.8 fb. 
The rate for e + e~ — > ttH° is always less than 1 fb for Ma > 120 GeV. Projected luminosities 
for a 1 TeV machine are around 200 fb -1 , which would give more than 400 tth° events/year. 

The dependence of the process e + e~ — > tth° on the energy of the Higgs boson is shown in 
Fig. |TT[ As tan j3 becomes large, the distribution is peaked at larger values of %h = 2E h / yfs. This 
is a kinematic effect due to the changing of the Higgs boson mass, since as tan (3 is increased, 
the Higgs mass increases with increasing Ma- 
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The production of a pseudoscalar in association with a heavy quark pair proceeds through 7 
and Z exchange, and also through the subprocesses, e + e~ — > A h® —> A°QQ, where h® = h°, H°. 
The rate is given by, M 
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Figure 8: Regions in the Ma— tan (3 plane where the cross section for e + e~ — > till®, (h® = h°, H°), 
production is larger than 0.75 fb at ^/s = 500 GeV. The upper left hand region results from 
e + e~ —>■ ttH°, while the region at the lower right is the result from e + e~ — > tth°. All NLO 
QCD corrections are included. The squarks are taken to have a common mass, Ms = 500 GeV 
and the mixing parameters are set to zero. 
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(16) 



where xa = 2Ea/^/s. Expressions for the coefficients, Fj(x^), are given in Appendix 0. 

The cross section for e + e~ — > ttA° is less than 10~ 2 fb at y/s = 500 GeV for all values of tan /3 
and Ma in the MSSM and so almost certainly will not be observable. However, the production 
of a scalar or of a pseudoscalar in association with a tt pair has been shown to be a useful 



mechanism for determining the CP coupling of the Higgs boson. [T3| In Fig. [12], we show the 
normalized distribution of the Higgs boson energy for the lightest scalar of the MSSM and for 
the pseudoscalar of the MSSM.f] The shape of the distributions is quite different, suggesting that 
this may be a useful method for discriminating between scalar and pseudoscalar Higgs- fermion 
couplings. For this to be practical, of course, requires a model where the pseudoscalar-fermion 
coupling is much larger than in the MSSM. 



5 These plots include contributions from the e + e — > h°A° resonance, which is specific to the MSSM. 
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Figure 9: Next-to-leading order result for e + e~ — > at yfs = 1 TeV using K = 0.94. The 
squarks are taken to have a common mass, Ms = 500 GeV. 

5 e + e~ — > 66/i° and e + e~ —> bbH° 

In the MSSM, the Higgs boson can also be produced in association with bb pairs and the 
differential cross section is given by Eq. @. The physics is significantly different from the 
production with a it pair, however. In the case of the b quark, the contribution of Fig. |3] 
receives a large resonant enhancement from the decay, A — > bb. This enhancement occurs 
when Mh. ~ Ma and Xh ~ 1 and so is relevant for Ma below about 120 GeV for e + e~ — > bbh°. 
We can estimate the resonant contribution in the narrow width approximation, 




a(e e -> N,h, ) | w ~ _ (— j (3— + — - 1 j , (17) 



where /3 is given by Eq.fl22|) in Appendix y and is the total decay width of the pseudoscalar 
A , shown in Fig. O. For the parameters considered here, the pseudoscalar is a very narrow 
resonance. This process is sensitive to the couplings of the Z to the A°h® pair in the MSSM, 6 



9zAh = cos(/5 - a) , 
6 The gzAhi couplings are normalized to the ZZhsM coupling. 
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Figure 10: Next-to- leading order result for e + e — > ttH° at y/s = 1 TeV using K = 0.94. The 
squarks are taken to have a common mass, Ms = 500 GeV. 



9zah = ~ sin(/3 - a) 



(18) 



along with the bbh® Yukawa couplings. 

Figs. O and [15] show the different contributions to the process e + e~ — ► bbh° for a low and 



a high value of tan/5, at y'i = 500 GeV [[ The curve labelled "iVW" is the narrow width 
approximation of Eq. (|T7D, labelled as "total" , while the curve labelled u AhZ n is only the 
contribution from the square of the diagram of Fig. At small Ma (< 120 GeV), the narrow 
width approximation is an excellent approximation to the total rate for these values of tan/3. 
For smaller tan (5, the narrow width approximation becomes increasingly inaccurate, since the 
Z exchange contribution becomes more and more relevant. Also at large Ma, the rate is given 
predominantly by the Z boson exchange contribution and is typically between 5 and 10 fb.0 
In Ref . , only the contribution from the Higgs coupling to the b quark, (the first diagram of 
Fig. [|), was included in the figures. This is only valid for Ma away from the resonance region. 
For Ma near the resonance region, Ref. || significantly underestimates the size of this process. 

7 We do not consider values of tan j3 smaller than 5 since the values of Mh corresponding to the resonance 
region of the bbh° production mode would be already experimentally excluded. 
8 In the Standard Model, the Z resonance is the dominant contribution. 
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Figure 11: Dependence of the process e + e~ — > tih° at yfs = 500 GeV on the energy of the 
Higgs boson, Xh = 2Eh/>/s, for Ma = 200 GeV. The squarks are assumed to have a common 
mass, Ms = 500 GeV and the scalar mixing parameters are set to zero. 



In the narrow width approximation, the QCD corrections to the rate are trivially included 
by including the QCD corrections to the pseudo-scalar width and this is done in Figs. [13] and 
|15| . The QCD 0(a s ) corrections to Ta can be found in Ref. 0. Away from the resonance, 
(large Ma), inclusion of the QCD corrections would require a complete calculation, which we 
do not include in the present analysis since the interesting region is near the resonance where 
the rate is enhanced. 

For the heavy Higgs production, the narrow width approximation is an excellent approxi- 
mation for all values of tan {3 so the QCD corrections can be accurately included everywhere. 
The rate for bbH° production is shown in Fig. 16. For tan (5 < 5, the cross section is larger than 
20 fb even for M A ~200 GeV. For tan/5>5, the rate is greater than 20 fb for Ma > 110 GeV. 
This process can potentially be used to probe the couplings of the heavier neutral Higgs boson 
and obtain a precise measurement of the product of the Higgs couplings, gbbHdzAH- 

The same type of behavior is still present if we vary over the possible values of the scalar 
masses and mixing parameters or if we look at y/s = 1 TeV. 
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Figure 12: Distribution in scalar energy for the pseudoscalar, A and the lightest neutral 
scalar, h° of the MSSM at = 500 GeV. The squarks are taken to have a common mass, 
M s = 500 GeV. 

6 e+e" bbA° 

We can safely work in the narrow width approximation also for the case of e + e _ — > 66 A 
production. In fact, in this case the h® resonances are completely dominant and the exact cross 
section can be distinguished from the one obtained using the narrow width approximation only 
at very high values of tan (3. The cases of tan (3 = 5 and tan (3 = 40 are illustrated in Figs. [17| 
and [TBI We can see that for tan (3 = 5 the exact cross section is indistinguishable from the 



one obtained in the narrow width approximation and the only relevant contribution is given by 
e+e" -> A ^ -> Abb. 

Unlike ttA° production, the process e + e _ — ► bb~A° is not suppressed relative to e + e _ — > 
66/i° production. The resonant behavior of the 66A° production mode is very stable with 
respect to the variation of the scalar mass and mixing parameters of the MMSM. Different 
parametrizations only affect the values of Mh and Mh at resonance. 
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Figure 13: Total width of the pseudoscalar, A , in the MSSM including NLO QCD corrections 
and assuming that there are no open channels for the A to decay into SUSY particles. The 
squarks are assumed to have a common mass, Ms = 500 GeV and the scalar mixing parameters 
are set to zero. 

7 Conclusion 

We have considered the production of the Higgs bosons of the MSSM in association with a 
heavy quark pair. The processes e + e~ — > tih°(H°) have rates around 1 fb at y/s = 500 GeV 
over much of the parameter space of the MSSM. They can reach rates of 2 fb in selected regions 
of the parameter space at y/s = 500 GeV and over most of the parameter space at y/s = 1 TeV. 
Observation of these processes could provide a precise measurement of the tth® couplings in 
a region not accessible at LEP2 and with 50 fb -1 , the tth° process could give a 90 % c.l. 
measurement with a statistical error of 

^ < ±10% (19) 

9tth 

over much of the Ma — tan (3 parameter space. The shape of the Higgs boson energy spectrum 
in these processes is quite sensitive to the value of tan (3. The rate for e + e~ — > it A is too small 
to be observed. 

For Ma > 120 GeV, the rate for e + e~ — > bbh° is also small, on the order of a few fb 
at y/s = 500 GeV. However, in the resonance region, Ma ~ M h , bbh° production receives a 
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Figure 14: Contributions to e + e~ — > bbh° at -^s = 500 GeV, for tan/3 = 5. The curve labelled 
'NW is the narrow width approximation of Eq. (|j~7D and it includes QCD corrections in the 
resonance region. The squarks are assumed to have a common mass, Ms = 500 GeV and the 
scalar mixing parameters are set to zero. 



large enhancement. This leads to the possibility of probing above the threshold accessed by 
LEP2. The rates for e + e~ — > bbH° and e + e~ — > bbA°aie greater than 20 fb throughout much 
of the MSSM parameter space and offer a unique window for studying the couplings of the 
pseudoscalar and heavier neutral Higgs bosons. 
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Figure 15: Contributions to e + e~ — > bbh° at y / ~s = 500 GeV, for tan/? = 40. The curve labelled 
'NW is the narrow width approximation of Eq. (|i~7l) and it includes QCD corrections in the 
resonance region. The squarks are assumed to have a common mass, Ms = 500 GeV and the 
scalar mixing parameters are set to zero. 

A Coefficients for e + e~ -> QQh°i Product ion 

The cross section for the process e + e~ — > QQh® (h® = h°, H°) is given in Eq. (|9]) in terms of 
coefficients, Gj. The coefficients G\ and G2 describe the radiation of the Higgs boson from the 
heavy quark (both photon and Z boson exchange) and are given by, @,[Jf] 



G 1 = 2 - ( - 4/3 [AMI - Ml) (2MJ + s) x h+ (20) 

(16M* + 2MJ - 2M 2 hi sx h + s 2 x 2 h - 4M 2 Q (3Mj - 2s - 2sx h )) A) , 



-2- 2 



Ga = 2 /^ QQ 2\ (- 96M Q + 2AM l M t ~ (~ M l + s + sx *) (-P 2 + x l)) + 

s 2 [p 2 — x 2 h ) x h \ 



x\) (24M* + 2 (Mj - W) + M% (-14M* + 12*r h + ^) ) A 



where we have defined 



9 Gi, ...G6 agree with those of Ref. B. G7 corrects a sign error in Ref. 0. We thank J. Kalinowski for a 
discussion of Ref. [fa . 
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Figure 16: Cross section for e + e~ — > bbH° at -^5 = 500 GeV, including QCD corrections in the 
narrow width approximation. The squarks are assumed to have a common mass, Ms = 500 GeV 
and the scalar mixing parameters are set to zero. 



A = log(^+f 

\Xh - P, 



Xt = [M 2 hi - M z z + s - sx h ) (21) 

and 

r f™2 ( rr min\2']\ rr max „ 1 n 1/2 

£={ { h U[X \ M , Xhl \ , (22) 

I „max „ j Q ) 

x h x h~r s 

with x% in = 2M h J^/s and x% ax = 1 - 4M^/s + M 2 Js. 

The other four coefficients, G3, . . . , G§ describe the emission of a Higgs boson from the Z-boson 
(plus relative interference terms) and can be written in the following form, 



G 3 = ]^ff^? (4i< + 12M| + 2M 2 z sx 2 h + s 2 (-1 + x h ) x\- 
Ml(8M z + s(-4 + 4x h + x 2 h ))) , 
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Figure 17: Total and partial contributions to the cross section for e + e~ — ► bbA at -^5 = 500 GeV 
and tan f3 — 5, including QCD corrections. The squarks are assumed to have a common mass, 
Ms = 500 GeV and the scalar mixing parameters are set to zero. 



(48M| + 12M£ - a (-24 + + 24x, - 1x\ 
4 9QQhi 9z M Q 



G 4 



Mz s * Xi S ( 6M| + Xh (~ M ^ ~ S + SXh ) ) + 

2 (Ml (M 2 hi - 3M| + a - sx h ) + M 2 Q (-4M 2 + 12M§ + 
-&9QQhi 9zM Q M z 



sxl)) A 



+ (4M< 



2 M 2 + 2s - sx ft ) A 



s 2 Xi 1 V <2 ^ 

where gr^ denotes the coupling of the Higgs boson to the Z boson 

gz = (V2G F ) 1/2 M z g ZZhi , 



where 



gzzh 
9zzh 



sin(P — a) 
cos(P — a) 



(23) 



(24) 



(25) 



The contribution from A h® production (see Fig. y) is given by 
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Figure 18: Total and partial contributions to the cross section for e + e~ — ■> bbA at -^5 = 500 GeV 
at tan (3 = 40, including QCD corrections. The squarks are assumed to have a common mass, 
Ms = 500 GeV and the scalar mixing parameters are set to zero. 



G 7 



9QQA9ZAhj / 2^(4^/2 _ 



Yi 



s 

+M 2 Q (4Ml-sxl))A 
where we have defined 

Y t = Ml -M 2 A + s-sx h 
while g QQA and g ZAhi are given in Eqs. (0) ([18]). 



sx h 



-9ZAhi9QQA- 



SX h — S — Mf 



hi 



Y 



2sx h (3(sx h - s - MQ -4,{M z h .{sx h -s-M l 



B Coefficients for e + e 

The coefficients of Eq. ( ]L6|) are0 



3 Our result for F4 (F3) is a factor of 2 (4) larger than the result of Ref. |§. 



(26) 



(27) 
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F i - 2 9qqa 
F 2 = 
+ 



x A - 



2M\ 

S 2 Xa 



A + 4 



2M 2 



M ^ ^ ' -24Mg + /3 2 S - S 4 

2 



sx A -Ml + 2Mq 
2 [ 2M\ + M 2 Q sx\ - 2M\ (Ml + sx A 



A 



s 2 /3 2 -x^ 
F 3 = 2/3 



(28) 



y 9QQh t gZAhj 



sxa — s 



M 2 A + AM 2 ^)(4.M 2 A - S x\ 



gQQAgQQhi9ZAh t 



- 2 [M\ + A*gsa& + Ml(-4M| + s - sx A )] A} , 



Q 



where we define 
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